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Abstract. Sand pile models are dynamical systems describing the evo- 
lution from N stacked grains to a stable configuration. It uses local rules 
to depict grain moves and iterate it until reaching a fixed configuration 
from which no rule can be applied. The main interest of sand piles re- 
lies in their Self Organized Cmticality (SOC), the property that a small 
perturbation — adding some sand grains — on a fixed configuration has 
uncontrolled consequences on the system, involving an arbitrary number 
of grain fall. Physicists L. Kadanoff et al inspire KSPM, a model present- 
ing a sharp SOC behavior, extending the well known Sand Pile Model. 
In KSPM(D), we start from a pile of N stacked grains and apply the 
rule: D— 1 grains can fall from column i onto the D— 1 adjacent columns 
to the right if the difference of height between columns i and i + 1 is 
greater or equal to D. This paper develops a formal background for the 
study of KSPM fixed points. This background, resumed in a finite state 
word transducer, is used to provide a plain formula for fixed points of 
KSPM(3). 

Keywords: Discrete dynamical system, self-organized criticality, sand 
pile model, transducer. 



1 Introduction 

Sand pile models were introduced in [T] as systems presenting a critical self- 
organized behavior, a property of dynamical systems having critical points as 
attractors. In the scope of sand piles, starting from an initial configuration of TV 
stacked grains the local evolution of particles is described by one or more iteration 
rules. Successive applications of such rules alter the configuration until it reaches 
an attractor, namely a stable state from which no rule can be applied. SOC 
property means those attractors are critical in the sense that a small perturbation 
— adding some more grains — involves an arbitrary deep reorganization of the 
system. Sand pile models were well studied in recent years (|8].[3].[5].|15|). 
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In [TT] , Kadanoff proposed a generalization of classical 
models closer to physical behavior of sand piles in which 
more than one grain can fall from a column during one it- > £) 
eration. Informally, Kadanoff sand pile model with param- 
eter D and N grains is a discrete dynamical system, which 
initial configuration is composed of N stacked grains, mov- 
ing in discrete space and time according to a transition Fig. 1: KSPM(D) 
rule : if the height difference between column i and i + 1 transition rule, 
is greater or equal to D, then D — 1 grains can fall from column i to the D — 1 
adjacent columns on the right (see figure [T]). 

In [10], the authors show that the set of reachable configurations endowed 
with the order induced by the successor relation has a lattice structure, in par- 
ticular it has a unique fixed point. 

More formally, sand pile models we consider are defined on the space of 
ultimately null decreasing integer sequences. Each integer represents a column 
of stacked sand grains and transition rules describe how grains can move from 
columns. Let h = (ho, hi, h®, ■ ■ ■ ) denote a configuration of the model, where 
each integer hi is the number of grains on column i. Configurations can also be 
given as height difference a — (co , o\ , o~i , . . . ) , where for all i > 0, Oi = hi — hi +x . 
We will use this latter representation throughout the paper, within the space of 
ultimately null non-negative integer sequences. 

Definition 1. The Kadanoff sand pile model with parameter D, KSPM(D), is 
defined by: 

— A set of configurations, consisting in ultimately null non-negative integer 
sequences. 

— A set of transition rules : we have a transition from a configuration a to a 

configuration a' on column i, and we note a A a' if 

• u\_ x = o-i-i +D - 1 (for i^O) 

• a[ = Oi- D, 

' a 'j = a 3 f° T 3 & ~ X ' * + D ~ 1 J'' 

Remark that according to the definition of the transition rules, a condition 
for a' to be a configuration is that <7j > D. We note a — » a 1 when there exists 
an integer i such that a A a'. The transitive closure of — > is denoted by A. A 
strategy is a sequence of nonnegative integers s — (si, . . . , st)- We say that a' 
is reached from a via s when a A a" A . . . A a' , and we note a A a'. We 
also say, for each integer t such that < t < T, that the column st is fired at 
time t in s (informally, the index of the sequence is interpreted as time) . For any 
strategy s and any nonnegative integer i, we state |s|, = #{t\st = i}. Let s , 

1 s° n s 1 i 

s be two strategies such that a —} er and a — > a . We have the equivalence: 

We say that a configuration a is stable, or a fixed point if no transition is 
possible from a. A basic property of the KSPM model is the diamond property. 
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If there exist two distinct integers i and j such that a -\ a' and a A a", 
then there exists a configuration a'" such that a' -4 a'" and a" A a'". From 
this property, one can easily check that, for each configuration a, there exists 
a unique stable configuration, denoted by 7r(cr), such that a —¥ tt(<t). Moreover, 
for any configuration a' such that a A a', we have the convergence property: 
7t(ct') = 7r(cr) (see [TU] for details). 

Let TV be a fixed integer. This paper is devoted to the structure of 7r((N, w )) 
(where W stands for an infinity of O's. n((N, 0^)) is denoted tt(N), for simpli- 
fication). Informally, our goal is to know what finally happens, starting from a 
configuration where all grains are in the first column. 

The main interest of our approach is to provide a new tool to study fixed 
points: a deterministic finite state word transducer, (a transducer is essentially 
a finite state automaton, which outputs a word during each transition). The 
idea is the following, we concentrate on a set I of D — 1 consecutive columns 
and constructs a set of states, according to the possible values of configurations 
on I. We then study how the way grains fall from the left part into / (input 
word) is related to the way grains fall from / to the right part (output word). 
The word transducer is formally defined in section [2~51 Using this transducer, an 
application to the case D = 3 is also proposed in section |2~51 This application 
leads to: 

Theorem 1. In KSPM(3), there exists a column number i(N) in 0(\ogN) such 
that: 

^(iV)[ iW ,co[ = (2,l)*[0](2,l)*0 w 
where o~uj] is the subsequence (<Xj, . . . , o~j), and [0] stands for at most one 0. 

This result can be interpreted as a kind of spatial emergence in a complex 
system. On a short length, the structure of the sand pile is complex, but a very 
regular shape is issued from the complexity. 

Describing and proving regularity properties, for models issued from basic 
dynamics is a present challenge for physicists, mathematicians, and computer 
scientists, there exists a lot of conjectures, issued from simulations, on discrete 
dynamical systems with simple local rules (sandpile model [3] or chip firing 
games, but also rotor router [12], the famous Langton ant [6] [7]...) but very few 
results have actually been proved. 

As regards KSPM(D), the prediction problem (namely, the problem of com- 
puting the fixed point) has been proven in [T5] to be in NC 3 for the one di- 
mensional case (the model of our purpose), which means that the time needed 
to compute the fixed point is in 0(log 3 N) where N is the number of grains, 
and P-complete when the dimension is > 3. In this paper we give a straight- 
forward characterization, describing asymptotically completely fixed points. A 
recent study ([S]) showed that in the two dimensional case the avalanche prob- 
lem (given a configuration a and a column i on which we add one grain, does 
it have an influence on index j?) is P-complete, which points out an inherently 
sequential behavior. 
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2 Avalanches and transducer 
2.1 Previous results about avalanches 

Let a be a configuration, o~^° is the configuration obtained by adding one grain on 
column 0. In other words, if a — (<7q, a\, . . . ), then o~^° = (ao+1, v\-> ■ ■ ■ )■ Clearly, 
for any strategy s, if <r A <r', then we have a"'' A ct'^ . In particular, since we 
have: (N,0 U ) A ir(N), we get: (A^O^ A ir(N)^ , i.e. (iV + 1,0") A 7r(A) i0 . 
The model is convergent, therefore we get the inductive formula: 

7r(7r(iV) i0 ) = ?r(iV + l). 

In the following, we will use the inductive approach described above, which 
consists in computing 7r(A + 1) by first computing n(N), then adding a grain 
on the origin column, and process all possible transitions until a fixed point is 
reached. For initialization, 7r(0) = W . 

The convergence property also allows to only consider leftmost strategies. A 
strategy s such that a A a' is called leftmost if it is the minimal strategy from 
a to a 1 according to lexicographic order. A leftmost strategy is such that at 
each iteration, the leftmost possible transition is performed. The k th avalanche 
s k is the leftmost strategy from 7r(fc — l)^ to ir(k). For our iterative approach, 
we need to describe avalanches. In a previous work |14j . we provide a simplified 
description of an avalanche. This description is a core result toward the construc- 
tion of the transducer in section 12.31 A first insight shows that any column is 
fired at most once within an avalanche. The formal statement of the avalanche 
description requires one more definition, even though it is graphically simple. 
For any sequence x = (x%, . . . , x n ) and any integers i, j with 1 < i < j < n, we 
denote xuji — (a^, . . . , Xj). 

Given an avalanche s k = (s\, . . . , s^.), a column s k is a peak if and only if 
s k > rnaxsj^ t _ 1 j. Intuitively, peaks are columns where an avalanche progresses 
rightward. 

Proposition 1. ml Let s = (si, . . . , st) be the k th avalanche and (pi . . . ,p q ) 
be its sequence of peaks. Assume that there exists a column I, such that for each 
column i with I < i < I + D — 1, i G ,s. Then for any column p such that 
p>l + D-l, 

p is a peak of s 7r(fc — l) p = D — 1 and 3i s.t. pi < p < pi + D — 1 

Furthermore, let pi — St, with Pi > I + D — 1, be a peak. Then 

T > t + pi — pi-\ — 1 and for all t s.t. t < t < t + pi — Pi-\ — 1, Sf = S{'— i — 1 

A graphical representation of this statement is given on figure [21 
The theorem and therefore the simplified description applies starting from 
a certain column / which depends on parameters of the model D and N. We 
say that the k th avalanche s k is dense starting at I and ending at m when m is 
the greatest fired column (V i > m, i ^ s k ) and any column between / and m 
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Fig. 2: Illustration of PropositionQ]with D = 6. Surrounded columns I to l+D — 2 
are supposed to be fired. Black column is the greatest peak strictly lower than 
£ + D — 1 before the avalanche. A column is grey if and only if its value is D — 1. 
Following arrows depicts the avalanche. 

included has been fired (V I < i < m, i E s k ). A consequence of Proposition Q] 
is that the avalanche s k considered is dense starting at I, where I denotes the 
parameter in the statement of Proposition!!] We define the global density column 
L(D, N) as the minimal column such that for any avalanche s k , with k < N, s k 
is dense starting at L(D, N). When parameters D and N are fixed, we sometimes 
simply denote L. 

Proposition 2. The global density column L(3,N) is in 0(\ogN). 

In KSPM(3), a trivial bounding of the maximal non empty column of ir(N) 
shows that it is in 0(^/N), so proposition Q] describes asymptotically completely 
avalanches used to construct the fixed point. We come back on this point in 
section [3j 

2.2 Successive avalanches 

When the k th avalanche is dense starting at I and ending at m, for each column i 
such that l + D — 1 < i < m, columns i — D + l,i and i+ 1 are fired within the k th 
avalanche. Therefore, ir{k)i = 7r(fc— — 7r(fc — 1)^. Moreover, n{k)j — Tr(k — l)j 
for j > m + D — 1. An intuitive consequence is that two consecutive avalanches 
are similar. This intuition is formally stated in this section. 

Let s k denote the k th avalanche of KSPM(D). We recall that the global 
density column L(D, N) of KSPM(Z?) is the minimal column such that for 
any avalanche s k , with k < N, s k is dense starting at L(D,N). We also de- 
fine ^(-D, N) — ((j) 1 , . . . , 4> n ), the subsequence of (s 1 , . . . , s ) reaching column 
L(D, N)+D — l. Formally, s k € <Z>(£>, N) -4=> L(D, N) + D - 1 e s k . ${D, N) 
is called the sequence of long avalanches up to N of KSPM(D). 

We also define the sequence (/jP, fi 1 , , /z ra ) of configurations such that /i° = 

7r(0) = W , and for each integer k such that <fi k = s m , we have /i k = 7r(m). 

The definition of long avalanche is motivated by the property above, which 
says that the effect of such an avalanche is easy to compute on the right of the 
global density column. 

Remark 1. In KSPM(D), if s k is a long avalanche, whose sequence of peaks is 
denoted by P k (the largest peak being max?'), from proposition Q] we have: 



- x(k) maxP k = n(k - l) maK pk -D+1 = 0, 
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- ir(k) t = vr(fc - 1), for L(D, N) + D - 1 < i < maxP fe , 

- 7r(fc)i = ir(k — l)i + 1 for maxP fe < i < maxP fe + D — 1, 

- Tr(k) l = n(k - l)i for i > ma,xP k + D - 1. 

This result is a clear application of transition rules (for each considered col- 
umn i we know which columns of the set {i — D + l,i,i + l} are fired in s k , so we 
can update the configuration). In other words, the main element that we need 
to compute (the right part of) 7r(/c) from 7r(/c — 1) is maxP fe . 

Lemma 1. In KSPM(D), let L be the global density column of N, and <P = 
(cj) 1 , . . . , 4> n ) its sequence of long avalanches up to N . Let k < n, and P k (resp. 
pk+i j fo e £fo e se q Uence of peaks i of cf> k (resp. 4> k+1 ) such that i> L + 2(D — 1). 
The largest peak of P k is denoted by max P k . We have: 

P fc \{maxP fe } = P fe+1 nI J L + 2(P-l),maxP' £ [ 

The lemma above can be seen as follows: |P fc+1 | > \P k \ — 1, and the \P k \ — 1 
first elements of P k+1 and P k are equal. Informally, the peak sequence can 
increase in arbitrary manner, but can decrease only peak after peak. 

Proof. Let k, k! be two integers such that (f> k is the K th avalanche, and <fi k+1 is 
the K lth avalanche. For each column i such that i G \L + D — 1, maxP fc [, we 
have i — D + 1, i, i + 1 e <j> k and therefore 7r(/i)i = tt(k — 

By definition of long avalanches, any avalanche s between <f> k and <fi k+1 stops 
before L + D — 1, i.e. for all i > L + D — 1, i ^ s. Combining it with previous 
remark, we have for all k" such that n < k" < k' 

for alii € lL + D-l,L + 2(D-l)l,n(n") l > tt(k - 1) (1) 
for all i e {L + 2{D - l),maxP fc [ , tt(k") 1 = tt(k - 1) (2) 

because columns within interval \L + D — 1,L + 2(D — I)[ can gain height 
difference when a column within [L, L + D — 1[ is fired. This is in particular true 
for k" = k' — 1. We now study the K ,th avalanche <j) k+1 . From relation ([1]) and 
since tt(k' — 1) is a fixed point, for alH € \L + D — 1, L + 2(D — 1)[, tt(k — 1); = 
D — 1 =>• tt(k' — l)i = D — 1. Let Q k (resp. Q k+1 ) be the sequence of peaks i of 
(j) k (resp. 4> k+1 ) such that L + D — 1 < i < L + 2{D — 1) Using proposition [1] we 
therefore get 

Q k C Q k+1 (3) 

Let I ={L + 2{D - 1), maxP fc [[. From relation © 

for all i G J, 7t(k - l)i = D - 1 it(k' - 1); = P - 1 (4) 

We now eventually prove the conclusion of the lemma. Let pi = min{i G P k }, 
from proposition [1] we equivalently have pi = min{i' G I\tt(k — l)j = D — 1} 
(the existence of is a hypothesis of the lemma). Let p' v = min-jY G P k+1 } = 
min{i' G I\tt(k' — l)i> = D — 1} (the existence of p\t is given by subset relation 
|3])), using relation Q we have p\, = p,-. 

Other peaks within / are obviously equal from proposition Q] and relation ^ 
with k" = «;'-!. □ 
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2.3 Transducer 



We now exploit the similarity between successive 
avalanches. Informally, we will cut configurations 
into intervals I\ , I2 , ■ • • of size D — 1 and study 
each of them and their interactions when consid- 
ering an avalanche. Given three successive inter- 
vals Ji-i, h and Jj+i, we construct a finite state 
word transducer which computes the influence 
of on knowing the influence of on 

Ii and the value of the configuration in Ii . The 
main idea to use transducers is that the value 
of any interval Ii with i > in 7r(0) is O 15-1 , 
so we can relate temporally emergent patterns 
arising from transduction iterations to spatially 
emergent patterns on stable configurations. 

The interval Ii is the column sequence ((D — 
l)i,(D - l)i + 1, ...,(£> — l)i + D — 2)). We 
call state of an interval Ii of a fixed point 

7T its value (TT(D-l)i,^(D-l)i+l--,^(D-l)i+D-2)- 

Hence, each interval state is an element of the 
set<S = {0,l,....,£)-l}' D - ;L . . 

We fix an interval Ii such that (D — l)i > 
L(N) + 3(D - 1). The largest peak j of (f> k , such 
that j < (D — is denoted by p(i, k). The type 
a(i,k) of the long avalanche 4> k on Ii is defined 
as follows. 

— if p(i,k) 6 /i-i, then a(i,k) = p(i,k) 
mod [D- 1]; 

— if fc) ^ ii-i, then a(i, fc) = e. 



in 



xD 



xc 



xm 



Therefore, the set of possible types is 
T = {e,Q,l,...,D - 2}. We say that two long 
avalanches are i- similar if they have the same 
type for i. Note that if a long avalanche <p k 

changes the state of Jj, then from remark [1] there m bold black. A light grey 
necessarily exists a peak of 4> k in the interval 
h-i. 

We now divide the sequence <P of long 
avalanches up to into maximal length subse- 
quences of the type (</> fe , fc+1 , 4> k ) such that, 
for each integer k' such that k < k' < k" , tfi k 
and 4> k +1 are i-similar. Such a subsequence is called an i-subsequence. An i- 
subsequence is said of type a for i when the type of each avalanche of the 
subsequence is a. When a is not the empty word e, we say that the subsequence 
is i-influent. Remark that, from Lemma Q] each (i + l)-influent subsequence is 



Fig. 3: D=4. Long avalanches 
up to 500, one by line. The 
global density column is lined 



square is a fired column, a 
dark grey square is a peak. 
The sequence of types of 4- 
influent subsequences up to 
500 is 0,1,2,0,1,2,0,2,1,0. 
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contained in an z-influent subsequence. See figure |3] for an example of i-influent 
subsequence. 

Lemma 2. Let $\k,k"] = {4^i <f> k+1 4> k ) be a subsequence of type a for i, 
with k" < n, and with Ii an interval whose columns are greater than L + 3(D— 1). 
Given the state (do, ax, cld-i) of Ii in the configuration /i fc-1 7 and a, one can 
compute, with no need of more knowledge: 

— the state (a' , a' 1; a' D _ 2 ) of Ii in the configuration fi k , 

— the sequence of types of the successive i + l-influent subsequences contained 
in $[k,k»] 

Proof. This is obvious when the type of the subsequence is e, since there is no 
change and the i + 1-subsequence contained in (4> k , . . . , (f> k ) is also e. 
The computation is simple when there is no integer m such that < m < a 
and a m = D — 1. In this case, the peak p(i,k) is the last peak of <fi k , thus 
A*p(j k) = 0' which gives that p(i, k) is not a peak of <f) k+1 , thus the subsequence 
is reduced to a singleton which is not i + l-influent (second part of the result). 
For (-D-l)i < j < p(i, k)+D-l, we have [i k = /i^ + l, and forp(i, k) + D-l < 
j < (D— l)(i + l), we have fi k = (J-j 1 - Thus, we have a' m — a m + l for < m < a 
and a' m = a m for a < m < D — 2 (first part of the result). 

Otherwise, (f> k contains a peak in Jj. Let q(i, k) denote the largest one. The 
column q(i, k) is the largest j such that fJtj 1 = D — 1 and j < (D — l)(i + 1). 
Thus q(i,k) mod D — 1 is the largest m such that a m = D — 1. In this case, <j) k 
starts an i + 1-subsequence of type q(i, k). Consider the following long avalanches. 
From Lemma [U while q(i, k) remains a peak of (f> k , p(i, k) also remains a peak 
of (f> k . From Remark [lj while q(i, k) is not the last peak of <fi k , the state of 
Jj remains invariant. So the first avalanche <p k that changes the state of Ii is 
the one whose last peak is q(i, k) (this avalanche exists from the hypothesis: 
k" < n). We have fi k , t fe j — 0, which closes the i + 1-subsequence of type q(i, k). 

We also have n k ' = /x* + 1 for q(i, k) < j < (D — l)(i + 1), and fi k — [i k for 

p(i,k) < j < q(i,k). This gives the state of Ii for [i k (as in the previous case, 
this can be rewritten to show that this state can be expressed only from a and 
(ao, ax, a.D-2)) an d proves that p(i, k) = p(i, k' + 1). 

The argument above can be repeated as long as we have a column j of Ii 
whose current value is D—l. When there is no more such column, the peak p(i, k) 
is deleted (its value becomes 0) by the next long avalanche which is necessarily 
<f) k from the maximality of t-similar subsequences. □ 

The algorithm below gives the exact computation. From the state of an in- 
terval Ii and an avalanche type on Ii, f returns the greatest fired peak in Ii, 
and g computes the new state of Ii and appends the result of / to a sequence 
of types on interval h+i. g recursively calls itself, anticipating the i-similarity of 
successive avalanches when maxP fc lies on the right of interval i. 
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" Input: a non empty type a and an interval state A = (ao, . . . , 00-2)- 
Data structure: a sequence w of types. 



Functions: 



f:SxT\{e}^T 


g : S x T \ {e} x T* -> S x T* 


f(A,a) := 


g(j4.,a,io) := 




if 


match /(A, q) with 




({m < a\a m =D- 1}^0) 


e -> (a + 1, . . . , a a + 1, a a +i, • • 


. ,ar>-2),w) 


then 


|p-> 




max{m \a m = D — 1} 


g((ao, • ■ • , Op-i, 0, a„+i + 1, . . . ,od 


_ 2 + 1), a,i«::p) 


else 

e 







. Computation: (A, a) \-¥ g(A,a,e) 



The algorithm above allows to define a deterministic finite state transducer 
C (see for example [2]). £ is a 5-tuple (Q, S, r, I, S) where the set of states Q is 
iS, the input and output alphabets (resp. £ and r) are equal to A = T \ {e} = 
{0, ...,£) — 2}, the transition function S has type Q x S — > Q X r* and is 
defined by the algorithm above: 5(q,a) = Computation(q,a). The initial state 
is (0, 0, 0), and we do not need to define a final state. The image of a word u 
by tT is denoted by t(u). 



If a is the type of 
an i-influent subsequence 
for a fixed integer i, then 
the sequence of types of 
the corresponding i + 1- 
influent subsequences ( i.e. 
subsequences where con- 
sidered avalanches lie) is 
t(a) . Thus, if u is the 
sequence of types of con- 
secutive i-influent subse- 
quences for a fixed integer 
i, then t(u) is the sequence 
of types of the correspond- 
ing i + 1-influent subse- 
quences. Note that the last 
considered avalanche may 
not be the last one of the 




last i + 1-influent subse- 
quence. Fig. 4: Transducer for D = 3 - Edges are labelled 
For the lowest interest- x\u, where x £ A is the type to the current inter- 
ing value, D — 3, the val (input) and u £ A* is the resulting sequence 
transducer £ can easily of types applied to the next interval (output). For 
be drawn. This diagram is example, t(abaaaaab) = abaab. Remark that, for 
given on figured]. For read- n > 0, we have: t((ab) n ) = (ab) 71 ^ 1 . 
ability, we write a and b in- 
stead of, respectively, and 1, for the alphabet of the transducer, and we omit 
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the drawing of states which are not connected with the initial one and are not 
useful for the computation of t(u), for any word u. 

The transducer has three transient states, (00, 10 and 20) and four recurrent 
states (11,12,21 and 22) organized in a cycle. A non trivial analysis of this 
transducer is given in appendix [X] The result is stated on the lemma below. 

Lemma 3. [D = 3/ For any k there exists n in 0(logk) such that for all u of 
length k, t n (u) is a prefix of (ab) u . 

2.4 From words to waves 

The lemma above can be used to describe fixed point configurations, noticing 
that | it | < N, as follows: 

Proposition 3. In KSPM(D). let L be the global density column of N and h 
be an interval whose columns are greater than L + 3(D — 1). Assume that the 
sequence of types of i -influent subsequences of long avalanches up to N is 



Let y be the size of the last subsequence of type p. We have y < x + 1, and 



Proof (sketch). It is a trivial induction on avalanches. We concentrate on the 
right part of fixed points: 7r(Ar)r^ j00 r. Initially for k = 0, it is equal to W . The D— 1 
first z-influent subsequences lead to D — 1, . . . , 1, W . Then from (D — 1, . . . , l) x , 
using lemma [T] to predict the size of each t-influent subsequence, we have that 
a sequence of type (0, ... , D — 2) corresponds to exactly (x + l)(D — 1) long 
avalanches, and the behavior verifies the following invariant : the ((x + l)p + y) th 
long avalanche, < p < D — 1 and < y < x + 1, has type p and lead to 



3 Conclusion 

Let us sum up results on KSPM(3). We introduced the finite state transducer 
which, given a sequence of types (associated with a sequence of long avalanches) 
on an interval J^, outputs the sequence of types (associated with the same se- 
quence of long avalanches) on interval h+i- We proved in a previous paper [14] 
that the global density column L(3,N) is in O(logiV), and therefore that to- 
ward the study of the fixed point ir(N), the word transducer applies starting 
from an interval Ij with j in J? (log N). Lemma |3] shows that iterating 0(log \u\) 
times the transducer function t outputs a prefix of (a6) w , from any input se- 
quence u. An upper bound for the size of any input word (sequence of types) 



(0,...,D- 2)^(0, . . . ,p), with x > and p < D - 2. 






□ 
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in KSPM(3) is N. As a consequence, there exists an index k in O(logiV) such 
that the sequence of types associated with subsequences of long avalanches up 
to N on interval Ij+k is a prefix of (ab) u . Finally, proposition [3] converts the 
temporal emergence of regularities when we iterate t into a spatial emergence 
of a wave pattern. It points out that as soon as a sequence of types which is a 
prefix of (ab) u is applied on an interval, then on the right of that interval n(N) 
is a wave. A simple framing of the maximal non-empty column of ir(N) shows 
that it is in 0(N), therefore the wave (2, 1)*[0](2, 1)* describes asymptotically 
completely fixed points of KSPM(3) obtained starting from a finite number of 
stacked grains. This concludes the proof of Theorem [T] 

We hope a generalization of this result to any parameter D, confirming ex- 
periments: 

Conjecture 1. For a fixed parameter D, there exists a column number i{N) 
in O (log AO such that: n(N) [i{N)>oo[ = (D - l,D - 2, . . . , 2, 1)*[0](D - 1,D- 
2,...,2,1)*0 W 

We name this pattern a wave for when your draw the corresponding sand 
pile, it looks like waves. Toward this aim, a possible outline is decomposed into 
two subproblems: one is to provide a general formula in 0(log N) for the global 
density column, allowing the use of transducers from that index; a second is a 
general study of KSPM(Z?) transducers resulting in the experimentally checked 
emergence of balanced outputs, then using proposition [3] we eventually conclude. 
Unfortunately from D = 4, transducers lack of human-readability for their num- 
ber of state is D 0-1 . Nevertheless, one may look for core properties on built 
transducers in order to deduce regular pattern emergence. 
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A Analysis of the transducer for D = 3. 

In this appendix we provide an analysis of the transducer for D = 3, leading to a 
proof of lemma[3] Note that though we consider maximal length subsequences of 
long avalanches, input words for the transducer may contain arbitrary numbers of 
successive occurrences of a and b since we consider only i-influent subsequences. 

We need some notations. Let q and q 1 be states of S and u be a word of 
A*. Consider, in the transductor, the path which starts in q, whose sequence of 
successive edge (left) labels is given by u. We say that we have qu = q' if this 
path terminates in q' . A word u is an entry word if 00 u is a recurrent state and 
for each prefix u' of u, 00 vl is a transient state. We denote by t q the transduction 
obtained changing the initial state into q. Hence ioo = t- We extensively use £21 j 
so we state ti\ — t' . A word u is basic for the state q if > 2 and for 

each prefix u' of u, \t q (u')\ < 2. For each current state q, the set of basic words 
for q and their images by t q are given below (trees represent case disjunctions 
according the beginning of u) 

(1, 1) : aaaa — » aba (2, 1) : aaa —> aba (1, 2) : aa — » ba (2, 2) : a — ► ba 
aaab aba aab — > aba ab — >■ ba b — > ba 

aab — > ab ab — > ab ba — > ba 

ab — > ab b —¥ ab bbu — > ab 

ba — > ba 
bb — > ba 

Each word u (such that t(u) 7^ e) admits a unique decomposition u = 
uqU\...u p such that uo is an entry word, for 1 < i < p, Ui is a basic word 
for the state 00 u$ui...Ui-i, and u p is a non-empty prefix of a basic word (for the 
state 00 UQU\...Up-i). The word u also admits a decomposition u — v! x u' 2 ...v! , 
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such that for 1 < i < p, Ui is a basic word for the state 21 u wi...Wi-i, an d u' p , 
is a non-empty prefix of a basic word (for the state 21 u' u' 1 ...u' pl 

A first result gives us a hint on the form of the sequence of types applied to 
successive intervals: 

Lemma 4. Let C be the language C = {abu, u G .4*} U {e, a}. 

— For each u G A* , we have t'(u) G C. 

— For each v G C, we have t(v) G C. 

— For each u G A* , we have t 2 (u) G C 

Proof. We prove the three items successively, using previous ones as hypothesis. 

— Let u G A* such that a/e. Consider the second decomposition seen above: 
u = u^u^.-.u'p,. We obtain t'(u) = t' '(u' 1 )t q (u' 2 ...u p ,), where q denotes a recur- 
rent state. 

• For p' > 2, t'(u[) is the image of a basic word for 21, thus t'(u[) G 
{ab,aba}, which gives t'{u) G C. 

• For p' = 1, t'(u) = t'(u[) and t'(u[) is the image of non empty prefix a 
basic word for 21, thus t'(w^) is a prefix of aba, which gives t'(u) G C. 

— Let i; G £. If v G {e, a}, then = e. Otherwise w can be written abu. Thus 

= t(ab)t'(u) = t'(u), and G C from the first item. This proves: 

t(v) G C. 

— Let u G A* such that e. We consider the first decomposition above: u = 
UqU\...u p . We obtain t(u) = t q (ui)t qui (u2...u p ), where q denotes a recurrent 
state. 

• For p = 0, £(m) = e, thus £ 2 (u) = e. 

• For p = 1, t(w) = t q {ui), and is the image by t q of a prefix of 
basic word for q, which gives that t(u) is a prefix of either aba or 6a 
(since possible images of basic words are ab,ba, and aba). This gives 
that t 2 (u) G {e, a}. 

• If p > 2, then t q (m) G {ab, ba, aba}. If t q {u\) G {a&, a&a}, then t(u) G £, 
thus t 2 (w) G £, from the second item. If t q {u\) = ba, then we can state 
t(u) = bau'. Thus t 2 {u) = t'{u'). We have t'{u) G C from the first item, 
thus t 2 (u) G £. 

□ 

Definition 2 (Height). T/ie height h of a finite word u G A* is h(u) = \\u\ a — 
\u\b\ where \u\ x is the number of occurrences of the letter x in u. 

Lemma 5. For any finite word v G C, we have: h{t{v)) < + 1 

Proof. This is obvious if v G {e,a}. Thus, stating v = abu, it remains to prove 
that, for any finite word u G A*, we have: h(t('u)) < + 1. 

Let us first consider the case when |m| — |u|f, > 0. Assume that we remove 
a pattern of the form ab or ba from u. This does not change the value of h(u). 
Moreover, for each recurrent state q, t q (ab) and t q {ba) both are elements of 
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{ab, ba} and qab = qba = q. This guarantees that pattern suppression does not 
change the value of h(t'(u)). 

Iterating this argument until there is no more pattern as above leads to the 
following fact: if we state u' = a h ^ u \ then we have h(t'(u')) = h(t'(u)). 

The integer h(u) can be written as h(u) — <U + r, with < r < 3. We 
have: t'(aaaa) — aba, and 21aaaa = 21. Thus t'(u') = (aba)* t' (a r ), which gives 
h(t'(u')) < h((aba) 1 ) + h{t(r)) <i + l. Thus h(t'(u)) < ^ + 1. 

The other case, when \u\ a — \u\b < 0, is similar. By simplifications of factors 
ba and ab, we obtain that h(t'(u')) = h(t'{u)), for u' = b h{u \ The value h(u) 
can be written as h(u) = 4j + s, with < s < 3. We have: t'(bbbb) = abbab and 
216666 = 21. Thus t'(u') = (abbab)n'(b s ), which gives h(f(u')) < h((abbab) j ) + 
h(t(s)) = j. Thus h{t'{u)) < + 1. □ 

Corollary 1. Given a word u € A* of length I, there exists an n(l) in O(logZ) 
such that t n ^(u) is a prefix of (a6) w . 

Proof. We first prove it restricting ourselves on words of C Given a finite word 
v on C, we define the maximal height g(v) = max{|/i(V)| \v' prefix of v}. The 
previous lemma gives the result g(t(v)) < 1 + ^r~- We can now use a trick to 
get the expected result. We define g'(v) = g(v) — |, then: 

9 (t(v))<i+ 9 M^mv))<^ 

From lemma 21 t(v) is element of C Thus we can iterate the inequality. By 
this way, we obtain, for each positive integer n: 

< ^ 

Thus, for n > \og 4 {g'(v)) - log 4 (|), we have: g'(t n (vj) < §, so g(t' n (v)) < 2 
and, by integrity, 

g(t n (v)) < i 

This last inequality enforces that u admits a decomposition t n (v) = Wi, w%. ..w q 
such that, for 1 < i < q, Ui € {ab, ba}, and w q £ {e,a,b}. Thus t n+1 (u) = 

t(wi)t' (W2), t'(w q ). Thus, t n+1 (u) is a prefix of the infinite word (ab)", since 

t'(ab) = t'(ba) = ab and t(ab) = t(ba) = e. 

Now, if we take a finite word u on A*, we have, from lemma 21 t 2 (u) € 
£. On the other hand, \t 2 (u)\ < 4|w|and \t 2 (u)\ + | > g'(t 2 (u)), which gives 
g'(t 2 (u)) < i\u\ + |. Therefore, for for n > log 4 (4|u| + |) — log 4 (|), we obtain 
that f n+1 (t 2 (u)) is a prefix of the infinite word (ab) u . In other words, for m > 
log 4 (4|u| + § ) - log 4 (|) + 3, t m (w) is a prefix of the infinite word (ab)" . □ 

Let us remark that I < N for any input word u so corollary [T] apply for actual 
sand pile behavior. We therefore have a strong property on words emerging from 
iterations of the transducer function t : they are exponentially quickly prefixes 
of (ubr. 



